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Finite volume scheme for two-phase flows in heterogeneous porous 
media involving capillary pressure discontinuities 
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Abstract 



o 
o 

04 ' We study a one-dimensional model for two-phase flows in heterogeneous media, in which the capillary pressure 

$_( ' functions can be discontinuous with respect to space. We first give a model, leading to a system of degenerated 

nonlinear parabolic equations spatially coupled by nonlinear transmission conditions. We approximate the 
solution of our problem thanks to a monotonous finite volume scheme. The convergence of the underlying 
discrete solution to a weak solution when the discretization step tends to is then proven. We also show, under 
assumptions on the initial data, a uniform estimate on the flux, which is then used during the uniqueness proof. 
A density argument allows us to relax the assumptions on the initial data and to extend the existence-uniqueness 
frame to a family of solution obtained as limit of approximations. A numerical example is then given to illustrate 

■ the behavior of the model. 

■ MSG subject classification. 35R05, 65M12 

, keywords, capillarity discontinuities, degenerate parabolic equation, finite volume scheme 

Introduction 

I The models of immiscible two-phase flows in porous media are widely used in petroleum engineering in order to 
^ ' predict the positions where oil could be collected. The discontinuities of the physical characteristics due to brutal 
I change of lithology play a crucial role in the phenomenon of oil trapping, preventing the light hydrocarbons from 

• reaching the surface. It seems that the discontinuities with respect to the space variable of a particular function, 
called the capillary pressure, are responsible of the phenomenon of oil-trapping [38, 10]. 

In this paper, we consider one-dimensional two-phase flows in heterogeneous porous media, which are made 
of several homogeneous submedia. A simplified model of two-phase flow within this rock is described in the first 

• section, leading to the definition of the weak solution. The transmission conditions at the interface between the 
] different submedia are written using the graph formalism introduced in [19] for the connection of the capillary 

pressures, which is simple to manipulate and allows to deal with any type of discontinuity of the domain, without 
any compatibility constraint, contrary to what occurs in [14] and to a lesser extent in [10, 23]. 

The graph way to connect the capillary pressures at the interfaces is well suited to be discretized by a monotonous 



H I Finite Volume scheme. A discretization is proposed in the second section of the paper. Adapting the material from 
the book of Eymard, Gallouet & Herbin [25] to our case, it is shown that the discrete solution provided by the scheme 
converges, up to a subsequence, to a weak solution as the step of the discretization tends to 0. The monotonicity 
of the transmission conditions is fundamental for proving the convergence of the scheme. 

Unfortunately, we are not able to show the uniqueness of the weak solution to the problem, because of the lack 
of regularity. As it will be shown in the fourth section, supposing that the fluxes are uniformly bounded with regard 
to space and time is sufficient to claim the uniqueness of the solution. The uniqueness proof is an adaptation of the 
one given in [19] to the case where the convection is not neglected. Here again, the monotonicity of the transmission 
conditions at the interfaces is strongly used. 

The existence of a bounded flux solution is the topic of Section 3. It is shown that if the initial data is regular 
enough to ensure that the initial flux is bounded with respect to space, then the flux will remain bounded with 
respect to space and to time. Such a result has already been obtained in [19], where a parabolic regularization 
of the problem had been introduced. A maximum principle on the flux follows. We also quote [10], in which a 
-BV^-estimate is shown on the flux. Since the monotonous schemes introduce some numerical diffusion, a strong 
analogy can be done between a uniformly parabolic regularization of the problem and the numerical approximation 
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via a monotonous scheme. The convergence of the discrete solution to a bounded flux sohition for regular enough 
initial data is thus naturally expected and stated in Theorem 3.1. The monotonicity of the transmission relations 
is essential during the proof. 

We are able to prove the uniqueness of the bounded-flux solution to the problem using the doubling variable 
technique. This work performed in Section 4 is summarized in Theorem 4.1 

In Section 5, a density argument allows to extend the existence and uniqueness frame to any initial data, using 
the notion of SOLA (Solution Obtained as Limit of Approximation). It is a more restrictive notion than the notion 
of weak solution, even if we are not able to prove the existence of a weak solution which is not a SOLA. The main 
result of the paper is given in Theorem 5.1, which claims that the whole sequence of discrete solutions built using 
the flnite volume scheme introduced in Section 2 converges towards the unique SOLA to the problem. 

Finally, a numerical example is given in Section 6. This example gives an evidence of the entrapment of a certain 
quantity of oil under the interface. 



1 Presentation of the problem 

We consider a one-dimensional heterogeneous porous medium, which is an apposition of homogeneous porous media, 
representing the different geological layers. The physical properties of the medium only depend on the rock type 
and are piece- wise constant. 

For the sake of simplicity, we only deal with two geological layers of same size. A generalization to an arbitrary 
finite number of geological layers would only lead to notation difficulties. In the sequel, we denote by = ( — 1, 1) 
the heterogeneous porous medium, and by fii = (—1,0), ^2 = (0,1) the two homogeneous layers. The interface 
between the layers is thus {x = 0}. T is a positive real value. 

We consider an incompressible and immiscible oil- water flow through fl. Writing the conservation of each phase, 
and using Darcy's law leads to: for all {x, i) £ fi^ x (0, T), 

(j>idtu ~ dx [lJ.o.i{u) {dxPo,i - Pog)] = 0, (1) 

- cj)idtu - dx [/iu;,i(w) {dxPws - Pwg)] = 0, (2) 

where (pi £ (0, 1) is the porosity of the porous media f2i, u is the oil-saturation (then (1 — u) is the water-saturation), 
fifi^i is the mobility of the phase P = w,o, where w stands for water, and o for oil. We denote by Pp^i the pressure 
of the phase /3, by pjs its density, and by g the gravity. 
Adding (1) and (2) shows that : 

dxq = 0, 

where 

q = -PwA'^) (dxPiv.i - Pwg) - Po.ii'^) {dxPo,i - Pog) (3) 

is the total flow-rate. For the sake of simplicity, we suppose that q does not depend on time, even if all the results 
presented below still hold for q G i?T^(0,T), as it is shown in [13, chapter 4]. 
Using (3) in (1) and (2) yields: 

(l)idtu + dx ( , ^"'''^"^ , . q + Xi{u){po - Pw)g - K{u)dx{Po,i - Pw,i)] = 0, (4) 

\PoA\u) + PwaW J 

where 

M'u-) = TT"^ r^- 

Po,i[u) + PwAu) 

One assumes that the capillary pressure {Po.i — Pw.i) depends only on the saturation and of the rock type. More 
precisely, {Poa — Pw,i) — T^i(u), where TTi{u) is supposed to be an increasing Lipschitz continuous function. The 
equation (4) becomes 

(j)idtu + dx {fi{u) - Xi{u)dxTTi{u)) ^ 0, (5) 

where 

ft{u) = — rTV~^~~rT9W + Xi{u){po - Pw)g- 

We do the following assumptions on the functions appearing in the equation. 
Assumptions 1 Fori — 1,2, one has: 
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1. TTi is an increasing Lipschitz continuous function; 

2. fioA is an increasing Lipschitz continuous function on [0,1], with /^o.i(O) = 0; 

3. fi^j is a decreasing Lipschitz continuous function on [0, 1], with /iuj.i(l) = 0. 

Remark 1.1 It is often supposed for such problems that the functions /i^ ^ are monotonous in a large sense, and 
that there exist irreducible saturations Si, Si £ (0, 1), with Si < Si, such that 

Aio,i(w) = i/ii e [0, Si], ^iw,i{u) ^Oifue [5*1,1]. 

// we assume that the functions ^p,i are strictly monotonous on their support, a convenient scaling would allow us 
to suppose that assumptions 1 are fulfilled. 

Wc denote by ^i{s) — \i{a)iT[{a)da, then (5) can be rewritten 

(jlidtU + dx (Mu) - da:(pi{u)) = 0. (6) 



Properties 1.1 It follows directly from assumptions 1 that for i = 1, 2 ; 

1. fi is Lipschitz continuous and fi{0) = 0, /i(l) = q; 

2. Xi is Lipschitz continuous, and Xi{0) = Xi{l) ~ 0, Xi{u) > ifu>0; 

3. Lpi is an increasing Lipschitz continuous fulfilling ^Pi{Q) = 0, ^p'i{Q) = '^[i'^) = 0. 

We deduce from the properties 1.1 that (6) is a degenerated nonhnear parabohc equation. 

Let us now focus on the transmission conditions through the interface {a; = 0}. We denote by a,; = hm^^^o i^iis) 
and I3i = hms-+i 7ri(s). We define the monotonous graphs tt; by: 

r 7r,(s) if.se (0,1), 

■^i('S) = N (-oo,ai] if s = 0, (7) 
[ [A,+oo) if.s = l. 

Let Ui denote the trace of u\q,^ on {x = 0} (which is supposed to exist for the moment). The trace on {x = 0} 
from Vii of the pressure Pjs.i of the phase (3 is stih denoted by P/3,i. As it is exposed in [23] (see also [19]), the 
pressure of the phase (3 can be discontinuous through the interface {x = 0} in the case where it is missing in the 
upstream side. This can be written 

mi("l)(^/3,l-^,3,2)+-/i,3,2(w2)(P/3,2-P;3,l)+ = 0, Pe{o,w}. (8) 

The conditions (8) have direct consequences on the connection of the capillary pressures through {x = 0}. 
Indeed, if < ui, U2 < 1, then the partial pressures Pq and Pw have both to be continuous, thus the connection of 
the capillary pressures 7ri(ui) = 7r2(u2) is satisfied. If ui = and < M2 < 1, then Pq^i > P0.2 and P^.i < Pw,2i 
thus 7r2(w2) < 7ri(0). The same way, ui = 1 and < it2 < 1 implies 7r2(M2) > 7ri(l). Checking that the definition 
of the graphs tti and tt2 implies 7ri(0) n 7r2(0) ^ 0, 7Ti(1) H 7r2(l) ^ 0, we can claim that (8) implies: 

Mui)r\n2{u2) ^9. (9) 

The conservation of each phase leads to the connection of the fluxes on {x = 0}. Denoting by Fi the flux in fii, i.e. 
for all X € ^li, 

Fi{x,t) = fi{u){x,t) - d^ipi{u){x,t), 
the connection of the fluxes through the interface can be written 

i^i(0,.) = F2(0,-), (10) 

where (10) has to be understood in a weak sense. 

We now turn to the problem of the boundary conditions. Because of technical difficulties occurring during 
section 4, we want that the solution to the flow admits bounded fluxes, at least for regular initial data. This will 
force us to consider specific boundary conditions, which will involve bounded fluxes. 
Let Gi ■ (a, b) 1— > Gi{a, b) (i = 1, 2) be a function fulfilhng the following properties : 
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• Gi is Lipschitz continuous, non-decreasing w.r.t. its first argument, and non-increasing w.r.t. the second. 

• for all a G [0, 1], Gi{a^ a) = fi{a). 

Let u, w e L°°(0,T), < u,u < 1 a.e., we choose the boundary condition 

= Gi(u(t),w(-l,t)), F2(l,t) = G2Kl,t),Ii(i)). (11) 

The way in which we approximate the boundary condition shall be judiciously compared with the discretization 
of the boundary conditions for scalar hyperbolic conservation laws using monotonous Finite Volume schemes (see 
[39]). 

We consider an initial data uq G L°°{fl), with < Uq < 1, then we can write the initial-boundary- value problem: 

(/jidtu + dj; lfi{u) - dj:ipi{u)] = in X (0, T), 
i^i(0,-) =i^2(0,-) on(0,T), 

^i(iii)n ^2(^2) 7^0 on(o,r), 

u(t = 0) = Uq in n, 

Fi{-l,t)=Giiuit),u{-l,t)) on (0,r), 
[ F2(l,<) = G2(w(l,t),u(i)) on(0,r). 



{V) 



We now define the notion of weak-solution 
Definition 1.1 A function u is said to be a weak solution to the problem (V) if it fulfills: 

1. ue L°^{nx {0,T)), with < u < 1; 

2. /or z= 1,2, ^,iu)eL\0,T■,HH^^)); 

3. for a.e. t £ (0,T), 7ri(ui(t)) n7r2(u2(t)) ^ 0, where ui denotes the trace of u^q. on {x = 0}; 
4-. for all ijj G I?(f2 X [0,T[), denoting by u{l, •) and u(— 1, ■) the traces of u on the boundary, 

/ / cj)iu{x,t)dt'>fi{x,t)dxdt + / cj)iUQ[x)^{x,G)dx 

/ [fi{u){x,t) - d^.Lpi{u){x,t)]dj;'il}{x,t)dxdt 

i{u[t),u{~l,t))^:{-l,t)dt- I G2(w(l,t),M(i))V'(l,t)di = 0. (12) 

Jo 

2 The finite volume scheme 

In this section, we build an implicit finite volume scheme in order to approximate a solution of {V) . We will adapt 
the convergence proofs stated in [21, 25, 23], which are based on monotonicity properties of the scheme. This will 
allow us to claim the convergence in L^iyi x (0,r)), up to a subsequence, of the discrete solutions built using the 
finite volume scheme towards a weak solution to the problem as step of the the discretization tends to 0. 

2.1 The finite volume approximation 

We first need to discretize all the data, so that we can define an approximate problem through the finite volume 
scheme. 

Discretization of il: for the sake of simplicity, we will only deal with uniform spatial discretizations. Let 
iV G N^, one defines: 

r Xj^j/N,^ VjGl-N,Nl, 
1 x,+„2 = ^-^^. VjGl-N,N-ll. 

One denotes by 5x = 1/iV. 

Discretization of (0,T): once again, we will only deal with uniform discretizations. Let M G N*, one defines: 
for all n G |0, M], = nT/M. One denotes by 5t = T/M. We denote by V the discretization of x (0, T) deduced 
of those oin and (0,r). 
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Discretization of uq: Vj E I-N,N - 1], 



Mo,i5(a;j+i/2) = "°+i/2 = g~ uo{x)dx. (13) 

Discretization of the boundary conditions: Vn G |0,M], 

= — / u{t)dt, = — / u{t)dt. 

dt (5t 7t„ 

The Finite Volume scheme: the first equation of (V) can be rewritten: 

(t)^^tu + d^F,{x, t) = 0, in x (0, T) 

with Fi{x, t) = fi{u) — dx^i{u). We consider the following implicit scheme: Vj e |— N, N — 1], Vrt G |0, M — 1], 

^.^!w/2_^fI±i/2^^ + FJ^ = (14) 

where F"^^ is an approximation of the mean flux through Xj on (t", t"^^), and i is chosen such that {xj, xj+i) C O^. 
This notation will hold all along the paper. We choose a monotonous discretization of the flux: Vj G |— N + 1, — 1] U 
|l,N-ll, VnG I0,M-11, 



where Gi is the same function as the one defined in (11). We also define 

= Gl(u"+\«_JV+l/2), Fji+' = G2(«Ar-l/2,Il"+^), (16) 



^0 - <-^n"-i/2'"o4 y^') 

2(o92«tM~092«t^)) 

where Uol"^,Wo2^ moreover satisfy 

^iKl')n7i2«,t^)7^0. (19) 



Remark 2.1 T/ie choice of the boundary conditions F!^'^^ has been done in order to ensure 

iF'l+^l < llGilloo < 00, \Fj:,+^\ < IIG2II00 < 00. 

Thanks to the following lemma, such a couple {uq\^,Uq^^) is unique in [0,1]^, thus the discrete transmission 
conditions system (17)-(18)-(19) is well posed. 

Lemma 2.1 For all {a,b) G [0, 1]^, there exists a unique couple {c,d) G [0, 1]^ such that: 

( 2(^i(c)-y>i(a)) 2{^2{h) - V2{d)) 

G,ia,c) -G2(d,6) j-^ , (20) 

[ 7ii(c)n^2(d) 7^0. 

Furthermore, (a, 6) c and (a, 5) ^ d are continuous and nondecreasing w.r.t. each one of their arguments. 
Proof 

For i = 1,2, 7r~^ are continuous non-decreasing functions, increasing on [7ri(0), 7ri(l)] and constant otherwise. Then 
we can build the continuous non-decreasing function A, defined by 

r K ^ M 

A • < 2 

' I P 1^ G2(7r2"^(p),6) - Gi(a,-n-~^(p)) -t- — ((^1 o^fj^i(p) - (pi(a) + (y92 o-S-2"^(p) - </52(fe)) ■ 
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For all p such that A(p) = 0, the couple (tt^ (p), tt^ {p)) is a solution to the discrete transmission conditions system 
(17)-(18)-(19). It is easy to check, using the monotonicity of the functions G,; that for all p < min7ri(0), A{p) < 0. 
Symmetrically, for all p > max7ri(l), A{p) > 0. Thus there exists such that A(p*) = 0. 

Suppose that there exists i such that G (7ri(0), 7ri(l)), then since ipi is increasing, A is increasing on a 
neighborhood of p^, and then the solution to the system (20) is unique. 

Suppose now that p^, ^ ljj(7ri(0), 7ri(l)). Either < mini7ri(0), then c = d = 0, or > maxi7ri(l), then 
c = d = 1, or p^, E [iTk{l),TTi{Q)] for k ^ I. We can suppose without any loss of generality that p^, G [7ri(l), 7r2(0)], 
then the unique solution to the system (20) is given by c = 1, d = 0. 

To conclude the proof of the lemma, it only remains to check that (a, 5) i-^ A is decreasing w.r.t. each one of its 
arguments, then the monotonicity of A and tt^^ ensures that (a, b) ^ c and (a, 5) i— > d are non-decreasing. □ 

2.2 Existence and uniqueness of the discrete solution 

We will now work on the implicit finite volume scheme given by (13)-(19) to show that this approximate problem 
is well-posed. 

Definition 2.1 Let N,M be two positive integers and T) be the associated discretization affix (0,T). One defines: 

^ _f zGL°°(r!. X (0,r)) / V(x,,a;,+i) Cl7„VnG I0,M-11, 1 
1 z\(^xj.xj+i)xit'-,t'-+^] is a constant J ' 

and 

Xv^{ze L^in X (0,r)) / 1, 2, Z|o.x(o,T) G . 

One defines uxi{x^t) G Xv, called discrete solution, given almost everywhere in (—1,1) x (0,T) by: for all 
JG I-N,N- 11, /or aHnG |0,M-11, 

J UTi{x,0) ^uo^v{x) = 'if{x,t) G {xj,Xj+i), 

I UT){x,t) = Uj'^i/a if {x,t) G {xj,Xj+i) x (i„,<„+i], 

where {u'j^ii2j o-f^ given by the scheme (14). 

The monotonicity of the flux i^"^^ w.r.t. ^11°'^^ us to rewrite the scheme (14) under the form 

H,+i/2 [u-Xl/,.u-^,„ {<Xln),^) = 0' (21) 

where -ffj+1/2 is continuous, increasing w.r.t. its first argument, and non-increasing w.r.t. all the others. 

Definition 2.2 A function vj) is said to be a discrete supersolution (resp. wx> is a discrete subsolution) if it 
belongs to X{T>), and if it satisfies: Vj G |— N, N — 1], 

(resp. i7j+i/2 (^ti/2^w!^+i/2 {<tl/2) < « )' 

Remark 2.2 A function u-p is a discrete solution to the scheme if and only if it is both a supersolution and a 
subsolution. 

Remark 2.3 It follows from the definition of the scheme, particularly from the definitions of the discrete boundary 
conditions (16) and of the discrete fluxes at the interface (17)-(18), that the constant function equal to 1 is a discrete 
supersolution, and that the constant function equal to is a discrete subsolution. 

We now focus on the existence and the uniqueness of the discrete solution to the scheme. In order to prove the 
existence of a discrete solution, we first need an a priori estimate on it. 
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Lemma 2.2 Let u-p be a discrete solution to the scheme associated to the initial data uq,©; iet v-d be a discrete 
supersolution associated to the initial data uo.-p, then for all t G [0,T], 

/ (j)t{uv{x,t) - vv{x,t))'^ dx < / (j)i{uo^Ti{x) - vo^-vix))^ dx. 

Symmetrically, if wx) is a subsolution associated to the initial Wo.t>, 

^ / (t)iiu-D{x,t) ~ 'w-D{x,t)y dx < ^ / (t)i{uo^x)ix) - wo,v{x)y dx. 

Proof 

Denoting by a~\b = max(a, 6), and a_Lfe = min(a, b), it follows from the monotonicity of the functions implies 
that 

^.+1/2 (^;:i/2-«}Vi/2T-;vv2 {<X\,2'^<Xln),^) < 0' 



where w-d is a subsolution. Since '"J_^i/2~'~^i+i/2 either equal to u^jXi/2 "-"^ '^'j+i/2 



^j + 1/2 I "j + 1/2 ' % + l/2'"i+l/2 I % + l/2 l^"fe+l/2 ' "'fe+1/2 



k^3 



< 0. (22) 



Thanks to the conservativity of the scheme, subtracting (21) to (22), and summing on j e |— N, N — 1] yields 

^ f ^^{uvix,t^+^)-Wvix,t^+^)y dx< I <l>^{uv{x,t'')-Wv{x,t^)ydx. 

Since this inequality holds for any n G |0, M — 1], it directly gives: Vi G [0, T], 

/ (l)i{uTi{x,t) - wt){x,t)y dx < y^ / (j)i{uT){x,0) ~ wv{x,0)y dx. (23) 



1=1,2'' ^'i i=l,2' 



The proof of the discrete comparison principle between a discrete solution and a discrete supersolution can be 
performed similarly. □ 
Let us now state the existence and the uniqueness of the discrete solution. 

Proposition 2.3 Let uq G L°°{Q), < uq < 1 a.e., then there exists a unique discrete solution uv to the scheme, 
which furthermore fulfills < up < 1 a.e.. Moreover, if stands for another initial data, Q < < \, approximated 
by Wo, I? following (13), and if we denote by v-p the corresponding discrete solution, then the following -contraction 
principle holds; 

/ (fii {uT)ix,t) - VT>{x,t))^ dx < ^ / (j)t{uo,T>{x) - vo,v{x))^ dx, Vte[0,T]. 

i=l,2'^^i 1=1,2"^^' 

Proof 

It follows from Remark 2.3 and from Lemma 2.2 that the following L°° a priori estimate holds: 

< uv{x, t) < 1, for all t G [0, T], for almost aU x G 51. 

Thanks to this estimate, mimicking the proof given in [24], we can claim the existence of a discrete solution u-p. 
Suppose that upi and vp are two solutions associated to the initial data uo,p and vo,-p. As it was stressed in the 
remark 2.2, both up and vpi are both discrete sub- and supersolutions. Then, Lemma 2.2 ensures that the following 
L^-contraction principle holds: 

Y / (l)iiup{x,t) - vp{x,t))''' dx < Y / (l3,{uo,p{x) - vo,pix))^ dx, VtG[0,T]. 

i=l,2'^^^> 1=1,2 -^^^i 

The uniqueness of the discrete solution up corresponding to the initial data Uq follows. □ 
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2.3 The L2((0,r);ifi(a)) estimates 

The current subsection is devoted to the proof of the discrete energy estimate stated in Proposition 2.4. Since 
the discrete solutions are only piecewise constant, we need to introduce discrete semi-norms, which are discrete 
analogues to the L'^((Q,T); H^{^li)) semi-norms. 

Definition 2.3 Let i = 1,2, one defines the discrete L^(0, T; semi-norms \ ■ ji.p.i on X'p,i by: Vz G X-p.i, 



^ll -D 4 = > ''^ > 



Sx 



n=0 jeJ,nt,i 

where J,nt,i = [-N + 1, -1] and J^a = P^N - l}. 

Proposition 2.4 For i = 1,2, one defines the Lipschitz continuous increasing functions 



■■ s ^ / Xi{a)TTlia)da. 
Jq 

There exists C > only depending on iti, (pi, T, Gi such that: 

E Mu'D)\lv,^<C. 
1=1,2 

This estimate is the discrete analogue to: 



E / / \d:^^^{u){x,t)\'^dxdt<C. 



1=1,2 

In order to prove Proposition 2.4, we will need the following technical lemma. To understand this lemma, first 
suppose that the total flow rate q is 0. Then, roughly speaking, it claims that, in the case where the capillary 
pressure is discontinuous at the interface, the discrete flux is oriented from the high capillary pressure to the low 
capillary pressure. Suppose now that q ^ 0. In order to respect the conservation of mass, some fluid will have to 
go through the interface, but we keep a control on the energy. 

Lemma 2.5 Let {a,b) G [0,1]^, and let {c,d) G [0,1]^ he the unique solution to the system (20), as stated in 
Lemma 2.1, then the following inequality holds: 

(-i(c) - ^2(d)) {G,{a, c) + ^^i^^^yf^) = (-i W - -2(d)) (c^id, b) + ^^^^^0^) > -klki(c) - 
Proof 

In this proof, we suppose that 7ri(0) > 7r2(0) and 7ri(l) > 7r2(l), the other cases do not bring any other difficulties. 
One has tti{c) n 7f2(d) 7^ 0, so there are three different cases: 

• 7ri(c) = TT2{d): in this case, one has directly: 

iMc)~Md))(G,ia,c) + ^l^^)=0. 

• (d) < 7i'i(0): the relation 7ri(c) fl TT2{d) ^ ensures that c = 0, thus it follows from the monotonicity of ipi 
and Gi that (^i(a) > (/?i(0) = 0, and Gi(a,0) > Gi(0,0) = /i(0) = 0. This gives: 

(-iW--2(rf))(Gi(a,c) + ^iM^)>0. 

7ri(c) > 772(1): this implies d = 1. From the monotonicity of ip2 and G2, we deduce that f2ib) < </'2(l) and 
G2(l,6) > G2(l, l) = q. This yields 

(^l(c) - 712 (rf)) (G2(d, b) + > - n2{d)\. 
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□ Proof of Proposition 2.4-- First check that the scheme (14) can be rewritten 



We multiply the previous equation by 5tTri{u''^^^^^^) and sum on j = —N, N — 1. This leads to 



where 



in+l 



N-1 
j=-N 



'] + l/2! l"i+l/2 "j'+l/2 



i^{-W,0,JV} 



C 



n+l _ 



1'i("J„l/2) (G** (""-1/2 '""1+1/2) ^ G'j(m"_^/2''""-1/2) 

1/2) (^'("j-1/2' "j+1/2) ^ ^'("^+1/2'""^+ 1/2) 



(5 a; 



n+l 



E ('^^k;i/2)-^^k-i/2), 

j<^{~N.O,N} 

StF^' (^i(^"i;2) - ^2«/V)) <5x - <5t^i(u'l|;2)/i(u^i;2) + (5t7r2«/V)/2«/+^) ; 
'5t^i(^i"1v'+i/2) (Gi (1 

+to2K^+_\/2) (G2 ("5 

Denoting by Lq a Lipschitz constant of both Gi 



E 



n+l 



U ,U 



'N + 1/2 



^1 I "-Ar+1/2' "-JV+1/2 



n+l -n+l 
*Af-l/2' " 



„ 1 ,,"+1 ,,"+1 

2 I "jV-l/2' "Af-1/2 



One has 



(7ri«|;,)-7r2K/V)) 



> -<5tiG(||7ri|U + ||^2||oo) 



+ ^^o"+'(^2K.l^)-^2K/V))- 



It has been proven in Lemma 2.5 that there exists Ci depending only on q and tt^ such that 

Using the definition of F^^^ , it is then easy to check that there exists G2 only depending on Gi, q, ni 

5x/2 

,,n+l 
-'0,2 



^"+1 > 5tC2 + 5t[ni{u'!X)2)~T^i{<f) 

+5t (7r2«/+^) - Tr2{uQ%^) 



'^2K/V)-^2«^^) 



5x/2 



Since tt^ is a non-decreasing function, : s i-^ I (f)iTTi(a)da is convex, then: Vn £ |0, M — 1] 



N-l 



^n+l> ^ (a.(u;;,%)-&.«+i/2))fe 



We denote by ^'j;(s) = / 'Ki{T)f[{T)dT, then an integration by parts leads to 
Jo 



- ^>,{a) = n,{a){G^{a,b) - f,{a)) - n,{b) {G,{a,b) - f,{b)) - / ttI{sM{s) - G,{a,b))ds. 



(24) 



(25) 



(26) 



(27) 
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Since /i(s) = Gi(s, s), it follows from the monotonicity of Gi and tt^ that 

b 

TT'^{s)(Ms)~G,{a,b))ds>0. 

J a 

Thus 



j^{-Af,0,Ar} 

So there exists C3, only depending on tt^, such that 

> (5tC3. (28) 

Let ii : s ^ i \J Xi{a)T:[{a)da, Cauchy-Schwarz inequality yields: V(a, 6) £ [0, 1]^, 
Jo 

(^.(a) - 7r,(6))(^,(a) - ^,(6)) > ^(a) - ^f- 

This ensures that 

V- (^^(";^i/2)-^^(";-i/2))' 

> d^t ^ ^ — — — ^ ; (29) 

Summing (24) on n G |0,M — 1], and taking into account (25), (27), (28), (29), (30), provides the existence of a 
quantity C, depending only on T, tt^, Gi, (pi such that 



(aw,r) -a(^n;2))' -6«t') 



16(^^)11.. + E ^ ^./2 ' + 5./2 I ^ (^1) 



=1,2 




□ 



Remark 2.4 T/ie estimate (31) is stronger than the one stated in Proposition 2.4, since it lets also appear some 
contributions coming from the interface. They will be useful in the sequel. Indeed, if we denote by u-p^i the trace of 
("f )|o interface {x ~ 0}, and if we denote by 7r>,j;(t) = Uq'^^ if t £ {n6t, (n + l)dt], then it follows from 

(31) that 

lim \\uv,t~-/v,t\\Lp(OT) =0, Vp G [1,00). 

Suppose that ux)^ converges in L^(0,T) towards a function Ui, as it will be proven later. Then, we directly obtain 
that jT>,i also converges towards Ui. Moreover, for all t > 0, 7ri(7_Dd,i(i)) ^ T^2{lDd;2{t)) 7^ 0- Since 

F = {(a, b) £ [0, 1]2 I ni{a) n n2ib) ^ (D} is a closed set of [0, l]^ 

we can claim that 7ri(ui) n7r2(M2) ^0 a.e. in (0, T). 



2.4 Compactness of a family of approximate solutions 

Let (Mp)pgN, {Np)p(zfi be two sequences of positive integers tending to +00. We denote T>p the discretization 
of X (0,T) associated to Mp, and Np. The L°°-estimate stated in Proposition 2.3 shows that there exists 
u G L°°{il X (0,r)). < M < 1, such that, up to a subsequence, u-p^, u in the L°°{yi x (0,T)) weak-* sense as 
p +00. 

We just need to prove that uxi^ u almost everywhere in x (0,r) to get the convergence of (wpp) towards 
u in L'"(0 x (0,r)) for any 1 < r < +00. To apply Kolmogorov criterion (see e.g. [12]) we need some estimates on 
the space and time translates of £,i{ux>)- 
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Lemma 2.6 (space and time translates estimates) For all 77 G M, for i = 1,2, one denotes fti^,j = {x G 
/ {x + rj) G i^i}, then the following estimate holds: 

muvK- + Vr) - C.iuv){-r)\\mn„,xio.T)) < Muv)\u-DMi\v\ + 2Sx). (32) 
One denotes Wi.v the function defined almost everywhere by: 

,„ I i^i^ry){x,t) m n,x{Q,T), 

w^,T>KX,i) - <^ ^ K2\(f}, X (0,T)). 

There exists Ci depending only on TTi, (j)i,T, d and C2 only depending on nt, (j)i,T, Xt, Gt such that: 

V77 e ||Wi,X'(- •)IU2(K2) < C177, (33) 

VTG(0,r), ||m,,2,(-,-+T)-m,,i,(-,0||L2(O.x(0,T-r)) <C2T. (34) 

The previous lemma is in fact a compilation of Lemmata 4.2, 4.3 and 4.6 of [25] adapted to our framework. The 
estimates (33) and (34) allows us to use the Kolmogorov compactness criterion on the sequence (wi^x>p)peN, and 
thus, there exists Wi G L?{VLi x (0, T)) such that for almost every {x,t) G (O^ x (0, T)), ^i{ux>p){x,t) — > Wi{x,t), 
and then thanks to the i°°-estimate < U'p^{x,t) < 1, one can claim that ^i(w-Pp) — > Wj in L'^{fli x (0,T)), for all 
r G [l,+oo[. Letting p tend to +00 in (32) insures that Wi belongs to L^{Q,T; H^{^i)). Since is a continuous 
function, we can identify the limit: 

Thus ^i(tt) G L^{0,T; H^{Qi)), and since ipi o is a Lipschitz function, there exists C depending only on 
T, TTi, (j)i, Gi, Xi such that: 

yt{u)\\mo,T;mi<n,)) < C, (35) 

and that S,i{ux>p) S,i{u), up to a subsequence, in L^iflt x (0, T)) a,s p +00 for any r G [1, +00). Since i — 1, 2, 
is an increasing function, one can claim that converges a.e. in fi x (0, T) towards u, and then: 

uVp ~* in the L°°{fl x (0, r))-weak- ★ sense, (36) 

uv, -> u in i'-(0, X (0,T)), Vr G [l,+oo[. (37) 

Roughly speaking, the approximation ud obtained via a monotonous finite volume scheme, which introduces 
numerical diffusion, is "close" to the approximation u'^ obtained by adding additional diffusion — eAw*^ to the 
problem. For such a continuous problem, we would have an estimate of type 

id^Uu'))dxdt<C\ 

which would lead to the relative compactness of the family (Ci(^*^))e>o in L^{fli x (0, T)). Then the family (^i(u'))^^Q 
is also relatively compact in i^((0, T); i?''(f2i)) for all s G (1/2, 1). This ensures that, up to a subsequence, the traces 
on the boundary and on the interface of (^i(u'^)) converge in L^{0,T). The continuity of and the i°°-cstimate 
ensure that the traces of u'^ on the boundary and on the interface converge in ^'^(O, T), for all r G [1, 00). 

This sketch has to be modified in order to deal with discrete solutions, which do not belong to i^((0, T); H''{Qi)) 
for s > 1/2. Nevertheless, a convenient estimate on the translates at the boundary, based on the discrete 
^^((0, r); iJ^(Oj;)) estimate stated in Proposition 2.4, leads to the following convergence result, which is proven 
in the multidimensional case in [17, Proposition 3.10]. 

Lemma 2.7 Let i = 1, 2, and let a G dfli. We denote by Uq,!?^,,; the trace of {uT>p)\ni '^^ ~ Q^}- Then, one has: 
for all r G [1, 00), 

Wa,-Dp,i "|Oi(a, •) in L'"(0, T) as p ^ +00. 

If we denote by Ui{t) = U|Q.(0,t), it follows from the remark 2.4 that 7ri(ui) n 7r2(u2) 7^ a-c. in (0,T). We can 
summarize all the results of this subsection in the following proposition : 

Proposition 2.8 Let (Mp)^, (Np)^ tend to 00 as p ^ 00, and let {T>p)^ be the corresponding sequence of discretiza- 
tions of fl X (0,T). Let (upp)^ be the sequence of corresponding discrete solutions to the scheme, then, up to a sub- 
sequence (still denoted by (^u-Op)^), there exists u G L°°{Vl x (0,r)), < u < 1 a.e., with S,i{u) G L'^{{{),T); H^{Vli)) 
(i = 1,2) such that: 

^ u a.e. in Vl x (0, T) as p ^ 00. 
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Moreover, keeping the notations of Lemma 2. 7, 



U-i^Vp.iit) u{~l,t) for a.e. t £ {0,T) as p ^ oo, 

Ui^x>p.2(t) u(l,t) for a.e. tG(0,T) as p oo, 

Uo.Vp.i{t) U|o^ (0,t) = Ui(t) for a.e. t e {0,T) as p ^ co, i = l,2, 

and 7ri(ui) n tt2{u2) ^ almost everywhere in (0,T). 

2.5 Convergence of the scheme 

Wc will now achieve the proof of the following result. 

Theorem 2.9 Let (Afj,)pgN, (iVp)pgN be two sequences of positive integers tending to + oOj and the asso- 

ciated sequence of discretizations of x (0,T). Then, up to a subsequence, the sequence (uu^) of the discrete 
solutions converges in L^iyi x (0,T)) for all r G [l,oo) to a weak solution to the problem {V) in the sense of 
Definition LI. 

Proof 

As it has been seen in Proposition 2.8, the discrete solution converges, up to a subsequence, towards a function 
u fulfilling all the regularity criteria to be a weak solution. In order to prove the convergence of the subsequence to 
a weak solution, it only remains to show that the weak formulation (12) is satisfied by the limit u. 

In order to simplify the proof of convergence of the scheme towards a weak solution, we will use a density result, 
which is a simple particular case of those stated in [22] . 

Lemma 2.10 Let a,5 G M, a < b, then: {tp G C;?°([a, 6])/V'' G C^{{a,b))} is dense m iyi'9(a,6), q G [l,+oo[. 

This lemma particularly allows us, thanks to a straightforward generalization, to restrict the set of test functions 
■0 for the weak formulation (12) to 

T = G v{il X [o,r))/9,0 G P((u=i,2a) X [o,r))}. 

Let ^Z" G T. For j G |— Np,Np — 1], n G |0,Mp — 1], we denote by V'j+1/2 ^ ^'(2;j+i/2: ^")- Assume that p is large 
enough to ensure: 

*-il2 = *"ri, Vne[0,Mp-ll, (38) 

One has also 

V-^^i/a^O, Vj G I-Np,Np-ll. (40) 

For j G I-Np, Np - 1], n G |0, Mp - 1], let us muhiply equation (14) by V'"+i/2'^*' J ^ I^^^P' ^^p ^ 

n G |0,Mp - 11, we get: 

J\/p-l Np-l A/p-1 Np-l 

E E ^^(";:i/2-";+i/2)^;+i/2'5a:+ E E (^'vv-^"+^)v^;-i-i/2-o, 

n=0 j=-Np n=0 j=-Np 

which can be rewritten thanks to (38), (39), (40): 

Mp-l Np-1 Np~l 

E E ^'<";:i/2(^';+i/2-'/';+%)^^- E </'^«°+i/2^^i/2^^ 

n=0 j=-Np 3 = -Np 

A/p-1 

+ E^* E Fr'i^u/2~ ^7^1/2) 

n=0 j<^{~Np,0,Np} 
A/p-1 A/p-1 

- ^ StF^2+'^i^l,f^) + ^ <5tF;j+V(l,t") = 0. (41) 

Ti=0 n=0 
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Using the definition of F"^^ , wc obtain 

Ap + Bp + Cp + Dp + Ep^ 0, (42) 
with, using the notation V'"jv_i/2 = V'"jv+i/2' ^^'^ ^7V+i/2 ^ "^N-i/i^ 

Mp-l Np-l 
n=0 j=-Np 

Bp = - '/''■"i+i/2'/'°+i/2^2; ; 

A/p-l 

n=0 j = -Np 

Mp-l Mp-l 
n=0 n=0 

Since /i^p : (a^,*) >— * st ^^^'^ (a^%^) ^ {^j^^j+i) ^ converges uniformly towards —dt4' as p — > cx), 

and since u-p converges in L^{n x (0, T)) towards u, 

\im Ap = -[ V / (j),u{x,t)dtip{x,t)dxdt. (43) 
Jo 

Thanks to the convergence in L^{fl) of 0) towards uq, we have 

lira Bp = - Y / '/'iWo(a;)'0(a;, 0)da;. (44) 

We have to rewrite Cp ~ Cp + Cp, with 

Mp-l 

= E E /':(";-i/2)(^i-i/2 - v^;+i/2), 

"=0 i^{-Wp,0,Afp} 

Mp-l 

= E ^ (Gi(«"+V2> ^;:i/2) - M^]-i/2)) (^7-1/2 - ^;+i/2)- 

Thanks to Proposition 2.8, the quantity Cp converges towards — / / fi{u){x,t)dxip{x,t)dxdt as p cxo. 

Concerning Cp, since 

Gi(uJ+j^/2' ""+1/2) ^ /»("j-l/2 '^^("j-1/2' "^'+1/2) ~ Gi(""-l/2' ""-1/2) ^ ^g\u"^1/2 ~ ""^-1/2 I ' 

^"-l/2-'0"+l/2 



ICp^l <iG||a:.'0||oo / / |5upp(x,i)|da:di, 



and since i~> ^ ^j, on (xj_i/2, 2:^+1/2) ^ ) is uniformly bounded by llSa,-!/;!! 

1^2 

where 5uT,p{x,t) = u"_^i/2 ^ "^-1/2 °^ (2^^-1/2,2:^+1/2) x (i",t"+i), (j e |-Np + l,Np - 1]), and otherwise. It is 
easy to check, thanks to the discrete L'^{{0,T); H^{ili)) estimates stated in Proposition 2.4, that Su-p tends to 
in L^{n X (0, T)) as p ^ 00. Then 

lim Cp = - Y fi{u){x,t)dxi^{x,t)dxdt. (45) 
P^°° Jo ~^2J^' 
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Since, using Proposition 2.8, ipi{u-Dp) tends to (pi(u) G L? {{Q,T) \ {^i)) in the L?{Vli x (0, T))-topology, one has: 
\im Dp = - f V / (p,{u){x,t)dl^ilj{x,t)dxdt ^ f f d.j;ipi{u){x,t)d^ip{x,t)dxdt. (46) 

The strong convergence of the traces, stated in Proposition 2.8 aUows us to claim that 

T 



lim = Gi{u{t),u{~l,t))ijj{~l,t)dt+ G2iu{l,t),u{t))iP{l,t)dt. (47) 

P-°° Jo Jo 

We can thus take the hmit for p -> cxj in (42), and it foUows from (43)-(44)-(45)-(46)-(47) that u fulfills the weak 
formulation (12). □ 

3 Uniform bound on the fluxes 

In this section, we show that, under some regularity assumptions on the initial data, there exists a solution with 
bounded fluxes. This existence result is the consequence of some additional estimates on the discrete solution, and 
will be necessary to get the uniqueness result of Theorem 4.1. 

Definition 3.1 ^ function u is said to be a bounded-flux solution to the problem (V) if: 

1. u is a weak solution to the problem (V ) in the sense of Definition 1.1; 

2. d^ipi{u) belongs to L°°{n^ x (0,r)). 

In order to get an existence result, we need more regularity on the initial data, as stated below. 
Assumptions 2 We assume that: 
1- 9.</?^(wo) e L°°{n,), < uo < 1; 

2. 7ri(uo,i) n 7i'2(uo,2) 7^ 0; where UQ,i is the trace of uq\^i. on {x = 0}, 

Theorem 3.1 Suppose that assumptions 2 are fulfilled. Let (Afp)pgN, {Np)p^-^ be two sequences of positive integers 
tending to +00. Let {uT>^)p^tq be the sequence of the associated discrete solutions obtained via the finite volume 
scheme (14), and let u be an adherence value of the sequence {uxip)p£H. Then u is a bounded flux solution to the 
problem (V) in the sense of Definition 3.1. This particularly ensures the existence of such a bounded-flux solution. 

All the section 3 will be devoted to the proof of Theorem 3.1. We only need to verify the second point in Defini- 
tion 3.1, because we have already proven in Theorem 2.9 that u is a weak solution. So the aim of this section is 
to get the uniform bound on the fluxes. Such an estimate can be found in [19] in the case where the convection is 
neglected. It is obtained using a thin regular transition layer between VLi and $72, and a regularization of the initial 
data uo. This technique was also used in [10] to get a _Bl/-estimate on the fluxes in the case of a non-bounded 
domain fJ, and for particular values of the data (which are supposed to be more regular). In this paper, we only 
deal with the discrete solution, which can be seen as a regularization of the solution to the continuous problem {V). 

We extend the definitions of the discrete internal fluxes (15)-(18) to the case n = —1, i.e. in the time i = 0. For 
all.7Gl-N + l,N-ll,jV0, 

F'-G(u" ) '^^K+i/2)-^^K-i/2) , . 

^3 - + • (48) 

Thanks to Lemma 2.1, there exists a unique couple (uq 1, Uq 2) solution to the system 

^0^,0 0^ "^iKl) -'/'l("°i/2) „,o ^ '^2(^i;/2) - '/52(?i[l,2) 

Fo" = Gi(u^i/2,4i) = G2(u[j,2,U?/2) '-j^ , (49) 

^i«i)n^2«2)^0- (50) 



Remark 3.1 Uq i and Uq 2 are given by Lemma 2.1, and so they are different of Uo,i and Uo,2- 
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Lemma 3.2 There exists C > depending only on uq, ipi, q such that 

max I 

jGl-^N+l,N-ll 



max_^ JF^\<C. 



Proof 

Since ipi{uo) is a Lipschitz continuous function, and ip^^ is continuous, wo|r2i ^ continuous function, and there 
exists yj+i/2 G i^j^^j+i) such that ^ *^o(?/j+i/2)- Then (48) can be rewritten 

F!/ = G,(wo(%-i/2),Mo(yi+i/2)) ■ 

Using the fact that dx-^iim) e gives directly: Vj £ |-N + 1,N - 1] \ {0}, 

\F°\ < max IIG^IIoo + 2max (||a:j,<y5,(uo)||L~(oo) ■ (51) 

The monotony of the transmission conditions tti{uq i) n7r2(Mg 2) 7^ ^ ^ind 7ri(wo,i) ri7r2(Mo,2) 7^ imphes that either 
1*0 1 > uo,i and Mq 2 > "0,2, or Uq ^ < uq^i and 2 ^ i''0,2- Assume for example that Uq i > uo,i and itg 2 ^ wo,2 
— the other case could be treated similarly — then one deduce from (49) that: 

/ ^^ V2 ("0(2/1/2)) - </'2 ("0,2) ^ ^0 / / / n ^ '/'l ("0,l) " '^1 ("O (^-1/2)) 

G2(mo,2, "0(2/1/2)) ^^jy^ < ^0 < Gi(uo(2/-i/2),wo,i) , 

and so since ipi{uo) is a Lipschitz continuous function, 

|Fo°| < max ||Gj||oo + 2max {\\dxfi{uo)\\L^(^n.)) . 

□ 

Proposition 3.3 There exists C > depending only on uq, ipi, d, such that 

max IF" I I < C. 

je[-N+l,N-ll Vne|0,M] ■' J 

Proof 

For aU j e |-N + 1, N - 1] \ {0}, for all n e |0, M - 1], 



] J 



{Gi{u'^j^ll2,u'^jXl/2) ^ ^i("j-l/2' ""+1/2)) + (^i("j'-l/2''""+l/2) ~ ("j- 1/2 ' "^ + 1/2)) 
/ (/5i(u":^i%) - <Pi(M"_i/2)\ / Vi(""+l/2) - Vi{u']+l/2)\ 



Thus, using (14) yields 



6x \ Sx 



5t ^i("j-l/2' "j + 1/2) ^'("i-1/2' "j + 1/2) 



pn+1 _ pn _ ut- "'V j-1/2' j + 1/2/ " ' v "j- 1/2 ' j + 1/2 / / j^n+l\ 



"j-1/2 "j-1/2 



r+h 

(j),5x -.."+1 ^"^+1 



(5t (""-1/2 '""+1/2) '^'(""-1/2' ""+1/2) ^j^„+i 



^^^■ + 1/2 "i+1/2 



y'»(M"+i/2) ~ '^i(^^"-i/2) 



n+1 Tpn+1\ 



"i-1/2 "j-1/2 
,n+l 



5t ^^{u"^l/2)-^^i^]-l/2) 



/pri+1 _ pn+l\ 



(jjiSx^ 



u 



n+l _ ri 
'i+1/2 "j + 1/2 



The monotonicity of the scheme is once again crucial, since it implies that there exist two non-negative values 
ct"j+i, a"j_^i such that 

(l + "w-i +Sj+i)^j ~ " Vi - ■ (5^) 
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The monotonicity of the graph transmission condition (19) ensures that either > Uq ^ and > ''^021 



4)^^ — "01 ^'^"^ ""o^^ — ""02- Suppose for example that Ugt^ — '"01 ^^"^ '^o^^ — ""02; the other case being 
completely symmetrical. 



^0 ~ ^0 



(Gi«|;2,<+i)-Gi(zi'!i/2,<i) 

G2K+\<;/)-G2K2,</2)) 
/ ^2Kl') - '^2K/2)\ / - <^2K/2)^ 



(5a;/2 



Sx/2 



(53) 



(54) 



It follows from (53) and from the monotony of Gi, ifi that 

- < (Gi«+;„<i) - G,{ul,„,ulS) + { 5^/2 ) • 

Similar computations as those done to obtain (52) provide the existence of a non-negative value a^'^J'i such that 

(1 + a-+_\) a-+\F-+^ < . (55) 

Considering (54) instead of (53) shows the existence of a non-negative value 69^^ such that 

(1 -f 6^+1) - b'^+^F^+^ > Fo". (56) 

We denote by j^^^ (resp. Jnihi^) ^^^^ integer such that 

Fl+\ = max (resp. F"+\ = min i^"+i). 



J„,.„ iel-N,: 



Either e {— iV, A^}, then it follows from the remark 2.1 that maXjg|_N_Nj < maxi=i^2 ||Gi||oo, or jjitx G 

|-N + 1,N - 1]. In the latter case, (52) and (55) imply 

maxF"+i = FZW < F\+i < max F". 



Similarly, (52) and (55) yield 



minF"+i = Fl+.\ > Fl+^ > minF" 



We obtain a kind of discrete maximum principle on the discrete fluxes, which corresponds to the uniform bound on 
the continuous fluxes proven in [19]. It follows from Lemma 3.2 that 



riG|0,Ml \jeI-N+l,N-l] 



< max ||Gi||oo + 2max (||a2;(/?i(wo)||L~(Oi)) 



□ Conclusion of proof of Theorem 3.1 Let (A^p)pgN, (Mp)pgN be two sequences of positive integers tending to 
-l-oo, and let {ux>j,)pen the sequences of associated discrete solutions. It has been seen in theorem 2.9 that (uT>p)p 
tends to a weak solution u in L'^{n x (0, T)), for all r S [1, +00). 

Let i = 1,2, let {x,y) £ fli, let t G (0,T]. For p large enough, there exists joi ji G "^int such that Xj„ < x < Xjg+i 
and Xj-^ <y< Xj-^+i, and there exists n such that t € (t",t""'"^]. 

\(piiu'D^){x,t) ~ (p,{uv^){y,t)\ = |v»(wJ„+\/2) - '^i(""i+\/2) 



< E |^«K':i/2)-^^^K-i/2) 
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Using the definition of the discrete flux (15): 



30 



I ^, {uv, )ix,t)-^, [uv, ){v,t)\< Sx G, {u^+l^, , 

We deduce from Proposition 3.3 that there exists C > 0, depending only on uq, ipi, Gi such that: 

30 

\'Pi{u-vJ{x,t) - ifi{uT)^){y,t)\ < ^ SxG <G{\x -y\ + 2Sx). 

Letting p tend towards +oo, i.e. Sx and 5t towards gives 

\^^{u){x,t)~^,{u){y,t)\<G\x^y\. (57) 
So we deduce from (57) that d^ip^{u) <= L°°{n, x (0,r)). □ 

4 Uniqueness of the bounded-flux solution 

This section is devoted to the proof of Theorem 4.1, which is an adaptation of [19, Theorem 5.1] to the case where 
the convection is taken into account. 

Theorem 4.1 If u,v are bounded-flux solutions in the sense of Definition 3.1 associated to the initial data uq, 
vo, then for all p G [l,+c)o[, u and v belong to G{[0,T]; LP(fl)), and the following -contraction principle holds: 

Vie[o,T], 

/ 4>i{u{x,t) — v{x,t))^ dx < / (f)i{uQ{x) — VQ{x))^dx. 
This particularly implies the uniqueness of the bounded flux solution to the problem (P) 

Obtaining a L^-contraction principle for a nonlinear parabolic equation is classical. We refer for example to [5, 26, 
34, 20, 31, 33, 11] for the case of homogeneous domains, and for boundary conditions of Dirichlet or Neumann type. 
We have to adapt the proof of the L^-contraction principle to our problem, and thus particularly to the boundary 
conditions and to the transmission conditions at the interface. 

We need to introduce the cut-off functions S C°'^(f2,]R+) defined by 

g , e - \x - a\ 

Lemma 4.2 For all £ P+([0,T)), 

liminf / 0{t) ^ / (^sign^{u - v){f,{u) - fi{v)) - {ipi{u) ~ ipi{v))^^ dxPo{x)dxdt > 0. 

Proof 

We define the subsets of (0,r) 

Eu>v = \ t e (Q,T) ui{t) > vi{t) or U2{t) > V2 

Eu<v = {Eu>vT = [te (0,T) I ui(t) < vi{t) and U2{t) < W2(i)} • 

Since the trace on {x — 0} of the function sign_|_(u — v){fi{u) — fi{v)) is equal to for all t G Eu<v, it is easy to 
check that 

lim / y2 sign^{u ~ v){ft{u) - fi{v))dxpldxdt = 0. 

Thanks to the fact that the trace of {(pi(u) — (pi(v))^ is equal to on the interface, one has also, 



liminf / 9 / dx {^i{u) — ipi{v))^ dxPodxdt > 0. 
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This particularly implies that 

liminf / 6* / sigii^{u - v){fi{u) - fi{v)) - {ipi{u) - ip,{v))'^ d^pldxdt>Q. (58) 

e^O /f - ^ — ' /o. L J 



Since u^v are two weak solutions, subtracting their corresponding weak formulation (12) for the test function 
tp{x,t) = 9{t)pQ{x) leads to 

/ X! / 4>iiu - v)pldtOdxdt + ^ / (j)i{uo - vo)pl9{0)dx 

+ [ [ - /«(^)) - C':.('P^(W) - '/'^(t'))) d^pldxdt = 0. 

Since tends to in L^(ri) as £ — > 0, one has: 

lim / ^ / (f), {u — v)pQdtOdxdt + ^ / (?!).j(ito - wo)Po^(0)d.i; = 0, 



thus 

lim / Oy^l ((/,;(«) - /<(«)) - d^(tp,(u) - (^j(w))) d^pgdxdt = 0. 
Thanks to the L°°{il x (0,T)) bound on the fluxes, one has 

^ X] / ~ /»(^)) ~ ^a;('^i(u) - d^p%dxdt 



(59) 



Li(0,T)i 



then, using a density argument, (59) holds for all 6 G L^(0,T). 

Replacing by Oxe^^^ in (59), and splitting the positive and the negative parts a = a"*" — a~, we obtain 

6* ^ / (sign_,.(M - v){fi{u) - fi{v)) - dx{fi{u) - fiiv))^) dxpldxdt 
^ / {&'^g'a--{u~v){f^{u)-f^{v))-dx{^i{u)-ifi{v)y)dxPldxdt +r(£). 



(60) 



with 



lim rie) = 0. 



For almost every t G Eu>v, it follows from the monotonicity of the graph relations for the capillary pressure 

7ri(ltl) n TT2{U2) / 0, ■^l(wi) n 7r2(w2) 7^ 0, 

that t belongs to Eu>v = |f G (0, T) ui{t) > vi{t) and U2(t) > W2(i)| ■ So we obtain exactly in the same way that 
for (58), that 

liminf/ ^ X! / sign_(u - w)(/,(u) -/;(«))- 9^ ((/?,(w) - V3i(u))~ dxpldxdt >0. 

It follows directly from (60) that 

liniinf / ^51/ sign_^(w - u)(/i(u) -/;(«))- 9^; (v3,;(ii) - V3j(v))"^ dxpldxdt >Q. (61) 

Adding (58) and (61) achieves the proof of Lemma 4.2. □ 
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Lemma 4.3 For all 9 e P+([0,T)), 



liniinf / 6 I (sign±(u - v) (Mu) - h{v)) - d^i^iiu) - M^))^) d^p^dxclt > 0, (62) 



/O Jf2i 

liminf / 9 [ (sigii±(u - v) (Mu) - /2(f)) - d,iip2{u) - ^2{v)t) d^p\dxdt > 0. (63) 



Proof 

For the sake of simplicity, we will only prove 



Jo. 



liminf j 9 j (sign+(M - v) (/2(u) - /2(f)) - dx{v2{u) - Lp2{v))^) dxp\dxdt > 0, 



but all the steps of the proof can be extended to the other cases. We denote by F^^y and Fu>v the subsets of (0, T) 
given by 

Fu>v ^{te (0, T) I uil, t) > t)} , ^,<„ - (F„>,)" ^{te (0, T) I uil, t) < v{l,t)} . 
Let e > 0. For almost every t G Fu<v, one has 

dx {f2{u){x,t) - ip2{v){x,t))^ dxp\{x)dx < 0. 

Then, using the fact that for almost every t e FuKv, the trace of sign_|_(u(-, t) — w(-,t)) {f2{u)(-,t) — /2(w)(-,t)) on 
{x = 1} is equal to 0, 

liminf / f (sign+(u - 1;) (^(u) - /2(f)) - dx{^2{u) - M^))'^) dxP\dx > 0. (64) 

We deduce from the weak formulation that for all 9 G 2?([0,T)), 

lim^ 9(^1^^ {f2{u)-dx^2{u))dxPldx-G2{u{l,t)Mt))^dt = 0. (65) 

Since the fluxes /2(u) — dxf2{u) and /2(w) — 9i:(/32(i') belong to L°°{^l2 x (0,2^)), a density argument, which has 
already been used during the proof of Lemma 4.2, allows us to claim that (65) still holds for any 9 G L^{0, T). So, 
it particularly holds if we replace 9 by 9xFu>^ This leads to 

lim/" 9!' [f2{u) - f2{v) - dx{ip2{u) - ip2{v)))dxp\dxdt = j 9{t){G2{u{\,t),u{t))-G2{v{\,t),u{t)))dt. 



It follows from the monotonicity of G2 that 

Vt G F„>,,, G2(u(l,t),u(i)) > G2{v{l,t),u{t)), 

thus 

liminf/ 9 [ {f2{u)- f2{v)^dx{<P2{u)-^2{v)))dxp\dxdt>Q. (66) 
In order to conclude the proof of Lemma 4.3, it only remains to check that 

liminf / 9 1 (/2(u) - f2{v) - dx{'f2{u) - ^p2{v)))dxp\dxdt 

(^sign_^(M - v){f2{u) - f2{v)) ~ dx{ip2{u) ~ ip2{v))'^^dxpldxdt. (67) 



lim inf / 9 



Since (^^^ is a continuous function, u{-,t) can be supposed to be continuous on Q2 for almost every t in (0,T). 
Particularly, for almost every t G Fu>v, there exists a neighborhood Vt of {x — 1} such that u{-,t) > v{-,t) for all 
a; G Vt. On Vt, one has 



R{x, t) = {j2{u) - f2{v) - dxi^2iu) - (P2(w))) - (sign+(M - v){f2{u) - f2{v)) - dx{v2{u) - ip2{v)y 
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= 0. 



Then, for almost every t G F^yv, 



lim 



R{x,t)d^p\{x)dx = 0. 



Moreover, since the fluxes f2{u) — dxf2{u) and f2{v) — dx'P2{v) belong to L°°{^2 x (0,T)), there exists C > not 
depending on e such that for almost every i, 



R{x, t)dxp\(x)dx 



0.2 

We deduce from the dominated convergence theorem that 



< c. 



lim 

e->0 



Jn, 



R{x,t)dxp\{x)dx = 0. 



This particularly implies that (67) holds. This achieves the proof of Lemma 4.3. □ 
Proof of the Theorem 4^.1. First, since u and v are weak solutions to a parabolic equation, they are also entropy 
solutions (see [26], [20]), and it has been proven in [18] that u and v belong to C([0, T], L^(ri)), in the sense that 
there exists u,v € C{[0, T],LP{n)) such that u ^ u, v — v almost everywhere in f2 x (0, T). 

Let u and v be two weak solutions, then some classical computations, based on the doubling variable technique 
applied on both the time and the space variable (see e.g. [26], [20]) yields yields that for any tp £ I?+(rii x [0, T)) 



(j)i{u{x, t) — v{x, t)) dtij{x, t)dxdt 

T r 
JUi 



(/)i(uo(a;) - vo{x))^'ijj{x,0)dx 



10 J^i 

sign^ (w(.T, t) - v{x, t)) {fi{u) {x, t) - fi{v) {x, t)) dxi^ix, t)dxdt 
/ dx{^Pi{u){x,t) - ipi{v){x,t))'^dxip{x,t)dxdt > 0. 

/O JUi 

Let 9 G P+([0, T)), then summing (68) with respect to i = 1, 2, choosing 

^(x,t) = eit) (1 - pii(x) - piix) - p\{x)) 

as test function, and letting e tend to leads to, thanks to Lemmata 4.2 and 4.3 : 



(68) 



\ dtQ(t)Y^ f Muix^t) ~v{x,t))'^dxdt+ f (f>,{uo{ 

^=1.2■^^^ 1=1.2-^^' 



x) -~vo{x))'^0{O)dx > 0. 



(69) 



Since u, v belong to C([0, T];L'^{n)), the relation (69) still holds for any € BV{0, T) with e{T+) = 0. Let t G [0, T], 
we choose 6 = X[o.t) in (69), obtaining this way the L^-contraction and comparison principle stated in the Theorem 
4.1. □ 



5 Solutions obtained as limit of approximations 

We aim in this section to extend the existence-uniqueness result obtained in Theorems 3.1 and 4.1 for any initial 
data uq G i°°(il), < uo < 1 a.e.. We are unfortunately not able to prove the uniqueness of the weak solution to 
the problem (P) in such a general case, but we are able to prove the existence and the uniqueness of the solution 
obtained as limit of approximation by bounded flux solution. Moreover, this limit is the weak solution obtained via 
the convergence of the implicit scheme (14) studied previously. 

Definition 5.1 A function u is said to be a SOLA (solution obtained as limit of approximation) to the problem (V) 
if it fulfils: 

• u is a weak solution to the problem (V ), 

• there exists a sequence {u^)i,i=n of bounded flux solutions such that 

u„ — > u m C([0, r]; i"'"(il)), asn^+co. 
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Theorem 5.1 Let uq G < < 1 a.e., then there exists a unique SOLA u to the problem (V) in the 

sense of Definition 5.1. 

Furthermore, if (Mp)pgN, {Np)p^jq are to sequences of positive integers tending to +oo, and if (M-Pp)pgN is the 
corresponding sequence of discrete solutions, then u in L^{fl x (0,T)), r G [1, +cxd). 

Proof 

The set 

f = {uo e I < lio < 1, d,ip,iuo) G L°-{n,), ^iKa) n ^2^,2) 7^ 0} 

is dense in {uq G | < uq < 1} for the L^(f7)-topology. Then wc can build a sequence (uq,^)^^^ such that 

hm ||uo.^ - woL WO) = 0. 

Let {uy)^ be the corresponding sequence of bounded flux solutions, then we deduce from the Theorem 4.1 that for 
all u, /.( G N, 

VtG[0,T], ^ / (j)i{u^{x,t) - Up,{x,t))^ dx < j (j)i{uQ^u{x) -Ui^,p{x))^ dx. (70) 

i=l,2^^> 1=1,2^' 

Then {uy)^ is a Cauchy sequence in C{[Q,T];L^{^)). thus it converges towards u G C([0, T]; L-'^(ri)). and 

ViG[0,T], ^ / (j)i{uu{x,t) - u{x,t))^ dx < j 4)i{uQ^y{x) - uq{x))^ dx. (71) 

1=1,2'' 1=1,2'' 

Let us now check that u is a weak solution. Since ipi is continuous, and since < < 1 a.e., Lpi{u^) 
converges in L^(r2i x (0,r)) towards (pi{u). The L^((0, T); 7/^(57^)) estimate (35) does not depend on uq, thus, 
up to a subsequence, {'Pi{u^))j^ converges weakly to (pi{u) in L^((0, T); 77^(17^)). It also converges strongly in 
L^{{0,T); H^{ili)) for all s G (0, 1). This particularly ensures the strong convergence of the traces of {ipi{u^))j^ on 
the interface. Since cp^^ is continuous, we obtain the strong convergence of the traces of (u^)^. Checking that the 
set 

F = {(a, 6) G [0,1]2 I ni{a)niT2ib) ^ 0} is closed in [0, l]^ 

the limits fulfill 7ri(Mi) n 7r2(w2) ^ 0, and so u is a weak solution, then it is a SOLA. 

If u and V are two SOLAs associated to the initial data uq and vq, we can easily prove, using the Theorem 4.1 
that 

VtG[0,r], j 4>^{u{x,t)-v{x,t)tdx<Y, j <l>^{.Uo{x)-Vo{x)fdx. (72) 

i=l,2'^^' z=l,2'^^» 

The uniqueness particularly follows. 

Let Uq G L°°(r2), < Mo < 1, and let {uq,^,)^ <Z E & sequence of approximate initial data tending to uq in L^(ri). 
We denote by u the unique SOLA associated to uq, and by [u^),, the bounded flux solutions associated to (uo.i/)^. 
Let (Mp)pgN, {Np)p^t>! be two sequences of positive integers tending to +00. Let p & v gN, let uUp the discrete 
solution corresponding to mq, and let u^^xip the discrete solution corresponding to uq.^. 

ll"I7p - '"llLi(f2x(0,T)) < ll^-D - U^,vJ\L^nx{Q,T)) + \\u,y.Vp - "i/ llii (Ox (0,T)) + Ikf " "II Li(nx (0,T)) • 

From the discrete L^-contraction principle (23), and from the continuous one (71), we have 

\Wvp - m||li(^^x(o,t)) < '^ min0- ^^""'^ ^ "o,i/,i7p ||li(0) + \Wi,.Vp - '"iv||Li(nx(o,T)) 

max</), 

+ T . Hmo.z. - M0||£i(O)- 
mm 0i 

Letting p tend to 00, it follows from the definition of (mq^i^.x)) (adapted from (13)) that 

lim \\uq x, - uo.,/,P„|Ui{o) = - uoIIli(o)- 

We have proven in the Theorem 3.1 that the sequence of discrete solutions converges, under assumption on the 
initial data to the unique bounded fiux solution, thus 

lim ||ui.,x)p - 'U'i^\\L\nxio,T)) = 0. 
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Figure 1: Saturation profiles for t = 20, t = 100, t = 200 



Tliis implies 
Letting v tend to oo provides 



II II max (p^ II II 

limsup \\uv^ - u||Li(nx(o,T)) < 2T . \\uo^^ - uoWl^q)- 
p—>oo mm (pi 



^lim^ \\uvp - '"|lLi(nx(o,T)) = 0. 



The convergence occurs in L^{0, x (0,T)), but the uniform bound on the sequence (u-Pp) in L°°(ri x (0,T)) ensures 
that the convergence also take place in all the x (0, T)), for p € [1, c»). □ 



6 Numerical Result 

In order to illustrate this model, we use a test case developed by Anthony Michel [32]. The porous medium fl ~ (0, 1) 
is made of sand for x G (0, 0.5) U (0.7, 1), with a layer of shale for x G (0.5, 0.7). 
First case: 

The total flow rate is equal to 0, since /sand(l) = /shaic(l) = 0, and the convection is the exclusive of the volume 
mass difference between the oil, which is lighter, and the water. The convection functions are given by: 

/sand(w) = 100 * /shale (w) = 50 * 



1 - 2m + 2m2 ■ 



The capillary pressures arc first given by 

7rsand(w) = , 7rshalo(") = 0.5 + 

The function v?sand and (pshaio, given by 

'Psand(u) = 10* / ^_^^^^^^2 ^Lndis)ds, (^shalc (w) = 0. 1 * 1 - 2 S + 25^ 



are computed using an approximate integration formula. The initial data uq is equal to 0, and u = 0.001, u = 0. 
The convection is approximated by a Godunov scheme, defined by 

min fi{s) if a < b, 
G,{a,b) = I «e[a,6] 

^ ' ' max ji(s) otherwise. 

se[b,a] 



A little quantity of oil enters the domain from the left boundary condition, and it moves forward in the first 
part made of sand. The discontinuity of the capillary pressure (figure 2) stops the migration of oil, which begins 
to collect at the left of the interface, as shown on the figure 1. One can check on the figure 3 that for t small 
enough, the oil-flux through the interface {x = 0.5} is equal to 0. The accumulation of oil at the left of {x = 0.5} 
implies an increase of the capillary pressure. As soon as the capillary pressure connects at {x = 0.5}, the oil can 
flow through the shale. The next discontinuity at {x = 0.7} does not impede the progression of the oil, since the 
capillary pressure force, oriented from the large pressure to the small pressure (here from the left to the right), 
works in the same direction that the buoyancy, which drives the migration of oil. 



22 



Figure 2: Capillary pressure profiles for t = 2Q,t — 100, t = 200 



Figure 3: Oil-flux profiles for t = 20, i = 100, t = 200 



For t ~ 200, the presented solution is a steady solution, with constant flux (figure 3). Some oil remains blocked 
in the first subdomain (0,0.5). Even if one puts u(t) = for t > 200, the main proportion of oil in the porous 
medium can not overpass the interface {x = 0.5} and leave the porous medium (0, 1). Indeed, the function defined 

by 

if a; ^ (0.4,0.5), 

'^saL(5(^-0.4)) if X £(0.4,0.5) 

is a steady solution to the problem for u = 0. It is easy to check that u{-, 200) > u^, thus the comparison principle 
stated in the Theorem 4.1 ensures that for aU t > 200, u{-,t) > ^i^ Thus for all t > 200 

^0.5 1-0.5 

u{x,t)dx > / u^{x)dx > 0. 



u%x) 



10 JO 

This quantity is said to be trapped by the geology change. Further illustrations, and a scheme comparison will be 
given in [32]. 

Second case: 

We only change the values of the capillary pressure functions (and also the linked functions v^sand and </'shaic). The 
amplitude of the variation of each function is reduced from 1 to 0.2, i.e. 

7rsand(u) = 0.2 * , 7rshalo(u) = 0.5 + 0.2 * u^. 

The graph transmission condition for the capillary pressure turns to 

(1 — '"sand)Wshalc = 0, 

where Ugand (resp Wshaic) denotes the trace of the oil saturation at the interfaces {x = 0.5} and {x = 0.7}. In this 
case, no oil can overpass the first interface, which is thus impermeable for oil. The only steady solution is 



u%x) 



1 if a; < 0.5, 
if a; > 0.5. 



An asymptotic study for capillary pressures tending to functions depending only of space, and not on the saturation 
has been performed in [13, Chapter 5&6] (see also [15, 16]). It has been proven that either the limit solution for the 
saturation is an entropy solution for th e hyperbolic scalar conservation law with discontinuous fluxes in the sense 
of [36, 37, 35, 1, 2, 3, 4, 6, 9, 8, 7, 27] (see also [28, 29, 30]), mainly when the capillary forces at the interface are 
oriented in the same direction that the gravity forces, or that non-classical shocks can occur at the interfaces when 
the capillary forces and the gravity are oriented in opposite directions. 
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Figure 4: Saturation profiles for t = 100, t = 500, t = 900 
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Figure 5: Capillary pressure profiles for t = 100, t = 500, t = 900 
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